Plug flow formation and growth in da Vinci Fluids 
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A new, da Vinci, fluid is described as a model for flow of dense granular matter. We postulate 
local properties of the fluid, which are generically different from ordinary fluids in that energy is 
dissipated by solid friction. We present the equation of flow of such a fluid and show that it gives 
rise to formation and growth of plug flow regions, which is characteristic of flow of granular matter. 
Simple explicit examples are presented to illustrate the evolution of plug flow regions. 
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On length scales much larger than grain size, dry sand 
appears to flow similarly to ordinary fluids. Apriori, it 
should be possible to construct continuum flow equa- 
tions for dense non-cohesive granular materials and in- 
deed such an approach has a long history. Yet, there 
is currently no agreement on any one set of such equa- 
tions as a clear favourite [ij . One of the main reasons for 
the complex behavior of granular materials is that the 
particles are inherently inelastic and energy dissipation 
dominates local dynamics. This feature makes conven- 
tional hydrodynamics ineffective and instabilities, such 
as inelastic collapse play a major role. As a result, 
large assemblies of macroscopic inelastic particles display 
a combination of both solid and fluid properties Q . The 
kinetic theory of so-called dense gases, which takes into 
considerations dissipation during collision between parti- 
cles, is only useful for low-concentration granular systems 
in very rapid flow 043- For dense flows, when many par- 
ticles rub against one another simultaneously, not only 
that formalism breaks down but also the concept of col- 
lision is not useful. Attempts have been made to extend 
the kinetic theory, introducing empirical stress tensors 
that take friction into account [SftUi], but the usefulness 
of this ad-hoc approach is not yet clear in the presence 
of multi-particle contacts. Much work is based on pro- 
posed phenomenological flow rules, relating strain rates 
and shear stresses [ij, [l^ i t>ut there is no agreement on 
any one such relation. Suggestions incorporating the idea 
of non-uniform stress propa gati on through force chains 
have also been put forward [14|, but the theor y b ehind 
these stress fields is not fully developed yet [15l - [l7| . 

A generic feature of granular flow is the formation and 
growth of plug reg ions, wherein the material moves as a 
rigid body Yet, most theoretical studies appear 

to assume existence of plug regions, rather than analyse 
the formation and evolution of such regions from flow 
equations. 



1rbbll@cam.ac.ukl Also at: Cavendish Laboratory, JJ Thomson 
Avenue, Cambridge CB3 OHE, UK 



The first purpose of this paper is to construct a minimal- 
model of flow for a dense system of non-cohesive grains 
from basic physical considerations, keeping to a simple 
mathematical description. Our second goal is to study 
quantitatively the formation and growth of plug regions 
directly from the equation. 

In this paper, we first derive an equation describing the 
flow of dense granular matter in regions of space, where 
the velocity field is not uniform. This involves the con- 
struction of a new stress tensor that captures the correct 
physics of particles interacting via hard core repulsion 
and solid friction. We then extend the model to include 
a description of the formation of plug regions and the 
dynamics of their expansion. The examples we present 
illustrate how local minima or maxima in an initial veloc- 
ity profile play the role of seeds around which plug zones 
form and expand. Our analysis shows that, at time t 
shortly after formation, the linear size of a plug zone is 
proportional to i^/'^. 

To construct our flow equation, let us consider first the 
difference between ordinary and granular fiuids. Concep- 
tually, the repulsive forces between grains are equivalent 
to interactions between molecules in ordinary fluids. The 
main difference is that in granular systems particles exert 
on each other non-central forces that cannot be viewed as 
stemming from two-body potentials and are due to solid 
friction [20|. Solid friction, described first by da Vinci 
[2l| . Amontons [l^ and Coulomb [1^, is an important 
energy-dissipating mechanism that transforms mechani- 
cal energy into heat stored in intra-granular degrees of 
freedom. This is different from the physics of ordinary 
viscous fluids. For example, the conventional dissipative 
term in the Navier-Stokes equation, which also represents 
conversion of mechanical energy to heat, only feeds en- 
ergy from macroscopic fluid disturbances to fluctuations 
on much smaller scales. Indeed, the viscosity term in the 
Navier-Stokes equation can be viewed as a result of drag 
forces on individual particles due to collisions with other 
particles, leading to viscous dissipation that is linear in 
the spatial derivatives of the velocity. The total force 
density, which gives the force on a volume element of 
the fluid due to other elements in its vicinity, is the sum 
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of two terms: a pressure gradient —VP and a dissipative 
term. In the following, we assume a simple dense granular 
fluid, which we call da Vinci fluid (dVF), after Leonardo 
da Vinci, who conducted the first recorded experiments 
on solid friction. The main difference between this and 
conventional fluids is that the dissipation is assumed to 
result only from internal solid friction, in contrast to the 
conventional drag forces that depend on relative velocity. 
The forces on a volume element in a granular fluid is the 
sum of the direct contact forces (normal and frictional) 
applied to its grains by grains in neighbor elements. We 
assume that the force density can be written as 

f =-V-a(")+Ff , (1) 

where the flrst term on the right hand side is due to 
normal intergranular forces applied on a volume element 
and the second is due to intergranular solid friction. It 
should be noted that the boundary between volume el- 
ements (which contain a very large number of grains) 
is a conceptual construct, and that the net normal in- 
tergranular forces between grains across a boundary is 
not necessarily normal to the boundary. Similarly, the 
net intergranular friction forces between grains across a 
boundary does not necessarily align tangentially to the 
boundary. 

The continuum solid friction term in a dVF can be de- 
rived in several ways. We have worked it out from 
detailed microscopic considerations, taking normal and 
solid friction intergranular forces into consideration. This 
derivation, however, is too long to reproduce here and it 
will be reported in detail elsewhere. Instead, we give in 
the following an alternative derivation from general ar- 
guments. Consider a region in the fluid, where the (non- 
symmetric) strain rate tensor Tij = diVj is non-zero (here 
and in the following, di and dt stand for d/dxi and d/dt, 
respectively). When the strain rate does not vanish, the 
friction force between adjacent volume elements should 
be proportional to the normal stress tensor, cr^"). The 
friction force should also be independent of the relative 
velocity between the elements ^3]. The simplest solid 
friction term consistent with this picture is — 7V-S', where 
S is the symmetric part of [cr'"' • T], 7 is proportional to 
an effective inter-granular solid friction coefficient /i, and 
T is a unit tensor, defined by Tij = Tij/ | T |, where | T \ 
is the norm of Tij . 

Let us consider first the equation of motion when the flow 
is nowhere spatially uniform. To make progress we write 
the 'normal' stress tensor in a simple form 

aJ)(r)^aS+P(p(r))<5ij , (2) 

where P is a scalar, p(r) is the local density and afj 
is position-independent, representing external and body 
forces. Since the material cannot be compressed to a 
point, the pressure has to diverge at some maximal den- 
sity pc- For the dense flows we aim to model, the average 



density, p, is not much lower than pc, 5p = p^ — p « Pc- 
A full description of the system is obained from: (i) an 
equation of state, P = P(p); (h) the equation of conti- 
nuity 

dtp + V ■ [pv] = (3) 
and (iii) from Newton's equation of motion, 

p [9tv + V Vv] = -V • (t("' - 7V • + f . (4) 

Here f is a local external force density acting on the sys- 
tem, e.g. due to stirring or interaction with the boundary. 
The flow equations ^ and (H]) evolve the density and ve- 
locity profiles of the system and they are valid wherever 
Vv does not vanish. In ordinary liquids, it is often pos- 
sible to replace the requirement of an equation of state 
by the simplifying assumption of incompressibility, but 
whether or not this is a good approximation for dense 
granular fluids is still an open question. 
The flow equations are complete and self-consistent, but 
they are still short of a full description of the behaviour 
of granular fluids. The reason is that under dense flow 
conditions non-uniform velocity fields are unstable due 
to formation of regions of uniform flow, called plug flow 
[3113: wherein Vv vanishes. Plug regions (PRs) are ex- 
pected to form and grow in the absence of stirring forces 
and our next goal is to obtain the flow equations for the 
motion of PRs, as well as for their formation and expan- 
sion. 

By definition, the acceleration field within a PR is spa- 
tially uniform; it is the total force on the region divided 
by its mass. However, it should be noted that changes in 
the position of the boundaries of the PR must be accom- 
panied by internal redistribution of stresses. The stress 
response rate does not enter the dynamics of the velocity 
field. We assume here that it is much faster than any 
other rate in the system and practically instantaneous. 
The total force on a plug occupying a region 17, due to 
the normal forces applied on the region boundary dU, by 
the adjacent fluid, is 

Fp = - / <7^"^(s) •n(s)ds , (5) 

where ds is an inflnitesimal surface element on 912 and 
n(s) is an outward pointing unit vector normal to the 

in) 

boundary at point s. cr^ is the stress just at the outer 
side of the boundary dO.. The friction force acting on the 
region is given by 

Ff. = -7/ S+(s) -11(8)018, (6) 

where S+ is the value of S also taken at the outer side of 
917. Due the distribution of intergranular contact orien- 
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tations at the boundary of a volume element, the bound- 
ary force density may have normal and tangential com- 
ponents, both of which are captured by S+. The total 
mass of the PR is 



Mp= [ p{r)d^r 
Jn 



(7) 



and, since the term v • Vv vanishes in il, the PR accel- 
eration is 



dtvp 



in 



/Mi 



(8) 



As we will demonstrate below, PRs are unstable and 
must grow. Therefore, to understand the flow behav- 
ior of dVF, we need to consider the kinematics of PR 
boundaries. To this end, it is convenient to define the 
scalar function 



(9) 



which is finite only outside PRs. The growth of a PR is 
equivalent to an expansion of the external contour lines 
of zero 'ijj surrounding it. The equation of motion for the 
scalar field ip is 



the initial conditions can change the density distribu- 
tion, which alleviates the need to solve for the equation 
of state. 

A straightforward example to analyse is when ^0(2;) is 
continuous and monotonic, i.e. dxVoi^) 7^ for all x. 
Then eq. ^ reduces to dtv = g, whose simple solution 
is that the velocity profile remains constant with time, 
v{x, t) = v{x, t = 0) + gt. 

A more interesting case is when the initial velocity profile 
has a local maximum, fixed at a; = (figure lb). A 
PR nucleates at the maximum, as described in detail in 
[23 ]. Basically, this is because the streamline at x = 
experiences friction forces opposing the fiow from both 
sides and it must decelerate relative to its surrounding 
fiuid. Eventually its velocity matches that of a neighbour 
streamline and they move together, forming a nucleus 
plug. The PR is slowed down by its surrounding fiuid and 
continues growing by the same mechanism. We wish to 
study the growth rate of the PR around the maximum. In 
regions not yet reached by the PR boundary, the velocity 
profile is monotonic and it does not change with time, as 
discussed above. Thus, for sufficiently short time t the 
growth of the PR is dominated by the velocity profile 
near the maximum. For convenience, we assume that 
the velocity can be expanded there as 



v{x,t — 0) = U — a{x/xo)'^ + O {x/xaf 



(12) 



The velocity profile at a later time t is given by 



where V^{y) is the velocity of the contour line of V'(r) at 

location r. Evidently, only the component of normal ""^^^ = Hx, t = O)+gt]0{t~ tp{x))+vp{t)0 (*p(^^ 

where 9 is the Heavyside step function, tp is the time 
when the plug boundary reaches point x and Vp{t) is the 
velocity of the PR at time t. The location of the PR 
boundary at time t is lp{t) (figure Ic) and the acceleration 
of the PR is 



to 9r2 is necessary to describe the expansion of the PR 
and this component is given by 



(11) 



The gradients on the right hand side of equation (fTTj) 
should be taken across the boundary and may be discon- 
tinuous. Nevertheless, a scrutiny of the numerator and 
the denominator will convince the reader that these dis- 
continuities and the corresponding 5-functions cancel out 
and leave a well-behaved term. This is also illustrated in 
the examples discussed in the following. 
To gain insight into the dynamics of growth of PRs in 
dVFs, we consider now an example that makes possible 
an explicit solution. Let a dVF be confined between two 
far-away boundaries at a; = ±L and between z = ±00 in 
the z-direction (figure la). We postulate an initial uni- 
form density p{x, t = 0) = po and an arbitrary initial 
velocity profile in the z-direction, vo{x)z — v{x,t = 0)z. 
Compressive forces in the x-direction are applied uni- 
formly across the boundaries at x = ±i. These forces 
give rise to a stress whose only non-zero component is 
axxj taken to be constant. Due to the symmetry of the 
system and the boundary loading, neither the forces nor 



ytVp 



Polpit) 



(14) 



The last term on the right of (ITi)) represents deceleration 
due to friction on the PR boundaries. The PR expands 
at a rate that can be found directly from (ITT]) . 



V* = dtlp = 



(15) 



Alternatively, p5)) can be obtained by solving for the 
time dt that it takes a streamline a distance dx away 
from the boundary to match velocity with the PR. From 
T5I) we obtain that the PR boundary grows as 



Ipit) 



2poVa 



1/3 



(16) 
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FIG. 1: a) A simple uni-directional flow in the z-direction un- 
der gravity g. The initial downward velocity profile increases 
from left to right, b) Same as in a but the velocity profile has a 
maximum, c) Nucleation and growth of a plug around a local 
maximum of the velocity profile. All the fluid accelerates at 
g except for the stream line through the maximum, which is 
slowed down from both sides by friction forces. Consequently, 
a plug forms and grows around this point. 

Eqs. (O, ([HI) and ^ provide a full solution for the 
velocity profile. The expansion of PRs as t^/"^, described 
by (fT6| . is generic and can be understood as follows. The 
friction force is proportional to the surface area of the PR 
boundary, while the mass is proportional to its volume. 
In our example, the friction force that the PR experi- 
ences from its surrounding material is constant, while 
its mass increases proportionally to lp{t) as it expands. 
Consequently, the acceleration decreases inversely pro- 
portional to lp(t). Similarly, given an initial velocity pro- 
file vq{x, y,t =: 0) that is cylindrically symmetric around 
a local maximum, the friction per unit height experienced 
by a PR of height h is proportional to 2Trlp{t), while its 
mass per unit height increases as ttZ^. The reduction 



in acceleration is then also inversely proportional to Zp, 
leading again to the result that Ip grows as t^^^. It can 
be shown that the linear size of the region (however de- 
fined) increases as t^^"^ for any velocity profile that can 
be expanded around a local maximum at (xq, i/q) as 

v{x -Xo,y-yo) -a[ 1 1 + ... . 

(17) 

To conclude, we have developed the flow equations of 
a fluid dominated by da Vinci - Amontons - Coulomb 
solid friction. The flow equations have been obtained for 
arbitrary flow regions, whether uniform or not. A key 
advantage of this model is that it gives rise naturally to 
unstable flows in the sense that the flow equations lead 
to formation and growth of plug regions. To our mind, 
this attribute of dense granular flow is extremely impor- 
tant but hitherto hardly studied. Our flow model is valid 
beyond these instabilities and it describes the formation, 
expansion and motion of the plugs region alongside the 
regions of nonuniform flow. We have discussed simple 
cases that are amenable to analytic treatment and we 
have shown that a generic feature of the flow is that, 
once a plug region has formed, its linear size increases 
with time t as t^^^. It is interesting that this growth law 
of PRs, which we have found in a model for dense granu- 
lar flows, has been observed also in simulations of muuch 
more dilute granular gases p5| . 

The phenomena described here resemble strongly obser- 
vations in flow of dense granular materials. Therefore, 
we propose this as a minimal model for such flows, when 
the material is dense and flows sufficiently slowly that 
grains maintain significant contact at all time. Experi- 
mental and numerical work, which we intend to take up 
in the future, is still needed to put this model to the test. 
In particular, it would be interesting to test the model's 
prediction for the growth rate of plug regions. 
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